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In this paper we consider nonlinear interaction between gravitational and electromagnetic waves 
in a strongly magnetized plasma. More specifically, we investigate the propagation of gravitational 
waves with the direction of propagation perpendicular to a background magnetic field, and the cou- 
pling to compressional Alfven waves. The gravitational waves are considered in the high frequency 
limit and the plasma is modelled by a multifluid description. We make a self-consistent, weakly 
nonlinear analysis of the Einstein-Maxwell system and derive a wave equation for the coupled grav- 
itational and electromagnetic wave modes. A WKB-approximation is then applied and as a result 
we obtain the nonlinear Schrodinger equation for the slowly varying wave amplitudes. The analysis 
is extended to 3D wave pulses, and we discuss the applications to radiation generated from pul- 
sar binary mergers. It turns out that the electromagnetic radiation from a binary merger should 
^et ■ experience a focusing effect, that in principle could be detected. 

<N 

>v 

S 

Recently much work has been devoted to the study of gravitational waves, largely due to the increased possibility 
of detection by facilities in operation such as LIGO (Laser Interferometer Gravitational Wave Observatory), or by 
ambitious detector projects under development such as LISA (Laser Interferometer Space Antenna) [lj . The possibility 
ly-j 1 of interaction between electromagnetic and gravitational fields has also led to alternative proposals for gravitational 
wave detectors, see e.g. Refs. 0-0 and references therein. Closer to the source, in an astrophysical context, the 
gravitational waves often propagate in a plasma medium, and the amplitudes can be much larger, which increases the 
number of possible interaction mechanisms, see e.g. Refs. |i||-|22|. Linear gravitational wave theory in a magnetized 
plasma has been studied by for example Refs including the backreaction from the plasma on the gravitational 

wave. In Refs. the authors have studied nonlinear responses to the gravitational wave by the plasma medium, 

O although the backreaction has been neglected. The nonlinear response gives raise to effects such as parametric 
instabilities 0, 0-IHHI3- large density fluctuations and photon acceleration 0. The application of gravitational 

wave processes to astrophysics has been discussed by for example Refs. -^H, and to cosmology by Refs. p^-[l8|. 
A number of works studying nonline ar prop agation of gravitational waves including the backreaction from the plasma 
have also been written, see e.g. \Wi IllL 119 ] . 

In Refs. pol l2l| geometrical nonlinearities from the Einstein tensor were considered, and a nonlinear evolution 
equation was derived. However, it was found that the nonlinear coefficient was proportional to the small difference of 
the phase velocity and the velocity of light in vacuum. In the present paper we neglect nonlinearities from the Einstein 
tensor and instead focus on the nonlinear response from matter. In particular we consider coupled gravitational and 
electromagnetic waves propagating in a magnetized plasma, with the direction of propagation perpendicular to the 
background magnetic field. For this wave coupling to be efficient, the interaction should be (almost) resonant, i.e. 
the electromagnetic wave propagation velocity in the plasma should be close to the speed of light. This increases the 
interaction strength and allows us to neglect the effect of the background curvature, in comparison with the direct 
interaction with matter. 

The plasma is modelled by a multifluid description, and we perform a self-consistent weakly nonlinear analysis 
of the Einstein-Maxwell system of equations for ID spatial variations in the high-frequency limit. This system is 
reduced to a single nonlinear wave equation for the coupled gravitational and electromagnetic waves. We then apply 
a WKB-approximation to this wave equation, and it turns out that the slowly varying wave amplitude obey the well- 
known nonlinear Schrodinger (NLS) equation |23|. In section HV Bl the analysis is expanded to include a 3D spatial 
dependence, allowing us to consider diffraction and/or nonlinear self-focusing of the wave. For certain conditions 
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the sign of the nonlinear coefficient is of focusing type, which for sufficient initial amplitudes implies solutions that 
undergo wave collapse. The conditions for collapse and the possible applications to astrophysics are discussed. 



The problem treated in this paper contains a number of small parameters used in the different approximations and 
expansions made throughout the paper. The aim of this section is to give an overview of these parameters, at what 
stages of the calculations they are introduced, and to what order the expansions need to be made. 

In section ITTT1 we introduce the small gravitational wave amplitudes, /i+,xi an d in calculating Maxwell and fluid 
equations, (|27|l - 1)33(1 . we only keep terms to (D(h+, x ). However, for the case when the gravitational wave interacts 
with an almost resonant Alfven wave, nonlinearities become important. The reason is that the resonance magnifies 
|<5B|/|Bo| where <5B and Bo are the wave magnetic field and background magnetic field, respectively. The magnification 
is of the order |<5B|/|Bo| ~ h + uj/5uj where 5ui is the frequency mismatch as compared to the Alfven mode. Therefore, 
we keep terms up to order 0(|<5B| 3 /|B O | 3 ), in the expansion of the wave equation (|41|) . in order to obtain the nonlinear 
amplitude modulation. Other nonlinearities originating from the Einstein tensor and the effective currents in l(7|l- l|12|) 
will not be magnified due to the resonance, and may be neglected in our treatment, see also the discussion following 



Furthermore, in section ITTll we introduce the high frequency approximation, valid for wavenumbers k 3> 1/fc, where 
r c is a characteristic radius of curvature of the background. This condition may also be written as the small parameter 
nB 2 /k 2 <C 1, where n — 8irG (we use units where c = fio = eo = !)• We have to keep terms to 0{nB^/k 2 ), in order 
to include the weak coupling between the gravitational mode and the Alfven mode. 

In section ITVl we first use the approximations <9 t « — d z , except in the linear wave operators, in order to simplify 
our system of equations. This is justified in the high frequency approximation. Furthermore, we use dt <C w c , i.e. we 
assume that the frequency of variations of fields is much smaller than the cyclotron frequency of all particle species 
in the plasma. In the expansion of Eq. (|36fl . we must include terms to 0(d t /uJc)- This is because the equations of 
motion to zeroth order in dt/u> c gives just the E x B drift of all particle species, and we must include terms of one 
order higher to obtain a nonzero current in the plasma. 

Next, in section TlVAl we introduce a weakly modulated wave (WKB-approximation), with the underlying assump- 
tion that the characteristic length scale for amplitude modulations obeys L|| 3> A. In deriving Eq. I(5()|l . we need terms 
to 0(X 2 /L 2 ) in order to include the dispersive term in (|50l) . We also use 1/C\ <C 1, where Ca is the Alfven velocity, 
and only keep terms to lowest non- vanishing order to simplify the equations. Then, in section TlV Bl we extend the 
problem to include small variations in the directions perpendicular to the direction of propagation. We apply the 
condition L± ^ A, where L± is the characteristic length for perpendicular variations. In subsequent calculations, we 
must keep terms to 0(X 2 *° obtain the diffractive term in Eq. I|67|l . see the discussion at the beginning of section 
II V Bl for more details. 

The parameters in this section are in principle independent of each other, with the sole exception that 8oj/uj ~ \jC 2 Al 
which is our main motivation for considering the regime C\ ^> I. In addition to these parameters there are a few 
other small parameters introduced throughout the paper, but those quantities are related in different ways to the ones 
presented in this section, and the relations are explained in the text. The only restrictions on the parameters in this 
section, besides their smallness, is due to physical applicability, i.e. the values of the parameters should be chosen to 
fit environments close to emitters of reasonably strong gravitational waves. An example of such parameter values is 
presented in the final section. 



II. OVERVIEW OF THE APPROXIMATION SCHEME 



Eq. (HOI. 



III. BASIC EQUATIONS 



The gravitational and electromagnetic fields are governed, respectively, by the Einstein field equations 



(1) 



and the Maxwell field equations 




= 3 
= 



;6 



(2) 
(3) 



V a Ff, c + S7 b F ca + V c F ao 



where G a b is the Einstein tensor, T a b is the energy-momentum tensor, F a b is the electromagnetic field tensor, j° is the 
4-current density and V denotes covariant differentiation. We use the spacelike signature ( — h ++) for the metric. 
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For the plasma we use a multifluid description, where the plasma is seen as a number of charged fluids, one 
for each plasma species, and collisions between particles are neglected. The total energy-momentum tensor is then 



Tab = T 



(fl) 



■ T 



(em) 



Here the fluid part is 
T (/0 _ 

1 ab " 



[(M(s) + P(s))u(s)aU(s)b +P(s)9ab] 



(4) 



where fA( s ),pr s ) an d u (s)a are t ne internal energy density, pressure and four- velocity for each plasma species, s, g a b is 
the metric tensor and the electromagnetic part is 



(em) 



ah 



= F a c F bc - y ab F cd F cd 



In the absence of collisions, the evolution equations for each fluid species can be written [2j 

Vrriab rriab ■ 
bT (s ) = F j b{s) 



(5) 



(6) 



as is consistent with Eq. 

With these preliminaries, we now follow the approach applied in Refs. 0, 0, 0, 0, |25|, |2(| and introduce an 
observer four- velocity, V a , so that the electromagnetic field can be decomposed relative to this into an electric and 
a magnetic part, E a = F a \,V h and B a — \t a bcF hc respectively. Here e a bc = V d e a bcd where t a bcd is the 4-dimensional 
volume element with £0123 = yjdet g a b\- Next we introduce an orthonormal frame (ONF) with basis {e a = e£d x i>}, 
where eo = V = V a e a , (i.e. V a = Sq), and write the fluid 4-velocity as u a — (7,7V), where 7 = (1 — v a v a )~^, 



a = 1, 2, 3, and v is the fluid three velocity. Dividing the four-current j a = 
Maxwell equations 0) and fluid equations © can be written 0, 0, |2^ | 



V-E = p + pe 
V-B = p B 
e E - V x B = -j - j B 
e B + V x E = -] B 

e o(7 n ) + V • {'ynv) = An 
(p +p)(e + v • V)7v = -7 



_1 Vp- 



■ 7v(e + v • V)p + qn(E + v x B) + (p + p)g 



in the same manner, the 



(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



where we have omitted the fluid species index, s, and it is understood that we have one pair of fluid equations <|lll 
1121 for each plasma fluid species. We have introduced the 3-vector notation E = (E a ) = (E 1 , E 2 , E 3 ) etc. and 
V = (ei,e2,e3). The charge density is p = J2s1(s)l n (s) an d n is the proper particle number density. The effective 
charges, currents and forces, originating from the inclusion of the gravitational field, are given by 



PE 
PB 

3E 

j-B 

An 



e a ^T%E 7 



-(r^-r| )^ 



T%E<*-e 



\ l f}0 "I 



E ^(r°£ 7 + rt£,) 



-7 [r£ + (r^ + r^ y 3 + r^V] 



(13) 
(14) 

(15) 
(16) 

(17) 
(18) 



where T^ c are the Ricci rotation coefficients associated with the tetrad {e a }. 

From now on we will assume that the plasma is cold, so that we can neglect the pressure terms in Eq. (|12JI and let 
u,( s \ = TO( s )ri( s ), where TO( S ) is the mass of each particle species. We will also apply the high frequency approximation 
[27| for the gravitational waves. Thus we will assume that the background gravitational field and the unperturbed 
plasma and electromagnetic fields fulfill Einstein's field equations (Q, and we introduce perturbations to the metric, 



9nv — 9$ +h fil/ , (greek indices are coordinate indices), and to the energy-momentum tensor, T a \, — T^' + ST a b, which 
should then fulfill 



(0) 



SG a b = nST ab 



(19) 
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where SG a b is the (linearized) perturbation to the Einstein tensor caused by the metric perturbation. 

It was recently shown by Ref. |l9j that in the high frequency approximation the gravitational wave can be taken to 
be in the transverse and traceless (TT) gauge even in the presence of matter. In this gauge the metric of a linearized 
gravitational wave propagating in the z-direction is given by |28| 

ds 2 = -dt 2 + (1 + h + )dx 2 + (1 - h + )dy 2 + 2h x dxdy + dz 2 (20) 

where h+ = h+(z,t) and h x = h x (z,t) denotes the two polarization modes of weak gravitational waves in the TT- 
gauge, and \h + \, \h x \ <C 1- The quantities h + and h x differ from the corresponding quantities in the vacuum case, 
h+ = h+(z — t),h x = h x (z — t), because of the weak interaction with the plasma. However, this difference is very 
small, which will allow us to use dt ~ —d z in many of the subsequent approximations. We note that components 
of the metric tensor corresponding to the quadratically nonlinear contribution from the pseudo-energy momentum 
tensor, as well as higher order cubic nonlincarities arc neglected in JSUJ.The basic motivation for this is that we 
neglect nonlinearities from the Einstein tensor, and focus on the nonlinear response from matter. When determining 
the regime of validity for this approach, a couple of things should be noted 

• As shown by Ref. |20j |. for uni-directional propagation of a weakly modulated plane wave in a flat background 
(in the absence of the gravitational wave), the nonlinearities in the Einstein tensor cancel up to cubic order as 
far as the evolution of the amplitude function is considered. (Although the background curvature due to the 
pseudo-energy momentum tensor of the wave is affected, this does not lead to amplitude modulations. We note 
that this result is consistent with the exact PP-wave solutions where the amplitude function is unaffected by 
the nonlinearities in the Einstein tensor, although the background curvature is modified due to the wave. ) 

• If the background in the absence of the gravitational wave is weakly curved (with a characteristic radius of 
curvature r c fulfilling kr c 1, where k is the gravitational wave number), nonlinearities in the Einstein tensor 
lead to amplitude modulations of the order of dth+, x ~ whS x / (k 2 r 2 ), see Ref. 21] . 

• If there are wave perturbations in matter close to resonance (i.e. assuming a wave mode in matter fulfilling 
lu = k + 5u>, where Su> <C u>), the matter response due to the gravitational wave is magnified. As will be 
shown below, for our case the nonlinear contribution, close to resonance, to the amplitude modulation scales 
as dth +iX ~ nToh^ x lu/5lu 2 , where To is the magnitude of the background energy-momentum tensor. We note 
that the nonlinear matter response dominates over that associated with the combined effect of the background 
curvature and the Einstein tensor nonlinearities, provided k~ 2 r~ 2 <C kTq/Suj 2 . 

• Using the high-frequency approximation (kr c » 1, s ee [27ll29| j1 and assuming a weakly nonlinear response from 
matter, it is always possible to use the TT-gauge As far as the direct interaction with the medium is 
concerned, i.e. the contribution that is magnified due to the resonance, the background metric can then be 
considered as flat (i.e. gfj = rj a b), see (2£j 

Next we choose our tetrad basis for the metric <|2(J|) as 

e = d t , ei = (1 - \h + )d x - \h x d y , e 2 = (1 + \h + )d v - \h x d x , e 3 = d z (21) 

Using H2(J|) and (|21|) in the linearized Einstein field equations (11911 . subtracting the n and 22, and adding the 12 and 
21 components of the equations one gets 0, El 113 

(d 2 -d 2 )h + = KiSTu-STn) (22) 
(d 2 -dl)h x = n{8T l2 + 5T 2l ) (23) 

In the basis (|21|l the non-zero terms in l|13l) - (|18l) are 

j E = -\{E 1 d t h + + B 2 d z h + + E 2 d t h x ~B 1 d z h x )e 1 

+ \{E 2 d t h + -B 1 d z h + -E l d t h x -B 2 d z h x )e 2 (24) 

j B = -±(Bidth + -E 2 d z h + + B 2 dth x +Eid z h x )ei 

+ ±(B 2 d t h + + E 1 d z h + -B 1 d t h x +E 2 d z h x )e 2 (25) 

g = -\i{v 1 d t h + +v 1 v 3 d z h + +v 2 d t h x +v 2 v 3 d z h x )e 1 
+ ^~f(v 2 d t h + +v 2 v 3 d z h + -vid t h x - viv 3 d z h x )e 2 

+ i 7 ((«? - v 2 2 )d z h+ + 2 Vl v 2 d z h x ) e 3 (26) 
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The generalization to a finite gravitational wave pulse width, i.e. including a weak spatial dependence perpendicular 
to the direction of propagation, leads to additional terms in (|24|) - (|26[1 and 122|) - (|23|) . An outline of this case is found 
in section llV Bl 

IV. COUPLED ALFVEN AND GRAVITATIONAL WAVES 

In Ref. |8j] a test fluid approach was taken to show that a gravitational wave can drive the amplitude of an 
electromagnetic wave to a nonlinear regime, and in Ref. [l9l | a linear analysis of the interaction between a gravitational 
wave and the extraordinary electromagnetic wave was made using the equations presented in section IIIII Here we 
intend to take a similar approach, but guided by the results in Ref. [8J, we will include a nonlinear response of matter 
and fields, while still assuming the metric perturbation to remain small. 

We assume the presence of a background magnetic field, Bn = -Bo e i, and introduce the perturbations: n = uq + Sn, 
B = (Bq + B x )ei, E = E y e2 + E z e^ and v = v y e2 + v z e% [3(|. We also note that in the case of gravitational waves 
propagating in a magnetized plasma, with the magnetic field perpendicular to the direction of propagation, only the 
/i+-polarization part of the gravitational wave couples effectively |3llto the electromagnetic wave, see also Ref. [l9j . 
Thus we put h x = in order to simplify the algebra. As in Ref. |8j, we assume v y <C 1, (because v y ~ h + ), and 
therefore neglect terms of the type v y and v y h + , but we allow for v z ~ 1. We will also consider slow variations such 
that dt uj c = qBo/m for each plasma species. With these restrictions the Maxwell and fluid equations ((7jl- ltT2*|l can 
be reduced to the following set of equations [32l |: 

d z E z = J2^( n o + Sn) (27) 

S 

d t E y -d z B x = -J2(qi(nn + Sn)v y )-^E y d t h + + ^(B + B x )d z h + (28) 

S 

d t E z = -Y,<ll{no + 5n)v z (29) 

S 

d t B x -d z E y = -\E y d z h + + \(B Q + B x )d t h+ (30) 

d t (7(^0 + Sn)) = -d z (7(710 + Sn)v z ) (31) 

d t ("/v y ) + v z d z ( 1Vy ) = ±(E y +v z (B + B x )) (32) 

d t ("fv z ) + v z d z (-/v z ) = ±(E Z - v y (B + B x )) (33) 

and the linearized Einstein equations, (|22(l - l|23|l . become 

{D 2 t - d 2 z )h + = n{E 2 y - 2B Q B X - Bl) (34) 

For the electromagnetic wave to exchange energy effectively with the gravitational wave, the two waves must be almost 
resonant, so we will consider all perturbations to be of the form B x B x (z — t), which allows us to use dt ~ —d z in 
simplifying our system of equations. A certain care must be taken here though; this approximation may of course not 
be used directly in the operator (df — d z ). Equations l|27|) - (|33|l can then be reduced to 

™ 2 ^p ( d z (jv 



(df - dl)B x + J2 ~^d z -fj-jf = d t {{B x + E y )d t h + ) + B d?h + (35) 
^ 9 \B + B X J 

9t(w>y)+v z d z (r/vy) = ±(v z (B + B x ) + E y ) (36) 



where uj v — y n$q 2 /m is the plasma frequency. In order to eliminate E y we can start with the linearized version of 
Eq. (|30|l . Using dt ~ — d z and integrating we obtain 

E y = —B x + \B h+ (37) 

Re-inserting this into (|30(l we find that the equation is fulfilled to 0(h+), which means that we can use the expression 
(|37|l for E y even when the electromagnetic amplitude is large, i.e. comparable to Bq. 
We now make an expansion of Eq. I|3()|l in the small parameter dt /u> c ancl use (|37|l to obtain v z . The result is 

Bq + B x ^ lo c ^ Bq + B. x ^ ^ 
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We see that v z does not depend on the charge of the particles, but is the same for all particle species, which means 
that the z-component of the current is zero, and the current is in the y-direction. Inserting v z , E y and expanding the 
wave operator in Eq. (|3"5f as d 2 — d z — (d t — d z ) (d t + d z ) ~ 2(d t + d z )d t we are left with the following coupled system 
of equations: 

{dt + V(B x )d z )B x = ^B d t h+ (39) 
(d 2 -d 2 )h + = -2kB B x (40) 

where V(B X ) = 1 - {1/2C%) (B / (B + 2B x )f /2 and we have introduced the Alfven velocity C A = (1/E, v%/uj 2 )V 2 . 
We note from that our previous assumption d t ~ —d z requires the background parameters to satisfy C\ ^> 1. 
In arriving at Eq. I|4UII we have also neglected a small term —nh+BoB x , resulting from the substitution of l|37l) into 
()34|1 . The left hand side of Eq. (|39|l contains the wave operator for compressional Alfven waves (also called fast 
magnetosonic waves), which in the considered regime propagates close to the speed of light. Furthermore, the left 
hand side of Eq. I|4(J|I describes gravitational waves in vacuum. The right hand sides of Eqs. (|39[l and l|4(J|) are 
the mutual interaction terms, which may provide a comparatively effective energy exchange, since the propagation 
velocities of the wave modes are close to each other. Now combining these two equations, and again expanding df — d z 
we obtain, after one time integration, the following wave equation for the combined electromagnetic and gravitational 
wave mode: 

(d t + d z )(d t + V{B x )d z )B x = -^B x (41) 



A. WKB-approximation for Quasimonochromatic waves 

Here we intend to show that for quasimonochromatic waves, the wave equation (|41|) leads to the nonlinear 
Schrodinger (NLS) equation [23 for the weakly varying amplitude. We begin the analysis of equation (|41|l by first 
considering the linear case. Thus replacing V(B X ) with V(Bo) = 1 — \/2C\ we get 

(dt + d z ) [d t + Jg.) 0.) B x = -^B x (42) 

From a plane wave ansatz, B x = Bexp(i(kz — u>t)), we then directly obtain the linear dispersion relation 

kB~ 

'A 



D(«,fc) = (u;-k)(w-fc+ 205-)- -^ = (43) 



This is in agreement with the dispersion relation presented by Ref. if we in their result make the appropriate 
approximations corresponding to d t ~ — d z . Solving the dispersion relation we get the following result: 



, ( 1 \ k 2 kB 2 

We see that there are two roots of l|43|) and we will hereafter refer to these two modes as the fast mode, for the root 
with positive sign, and the slow mode, for the root with negative sign. We see that for large k, i.e. k ^> yj kB 2 C a , 
we can neglect the small quantity kB 2 /2, and the dispersion relation for the fast mode takes the form u) » k, while 
for the slow mode we get u ~ k(l — 1/2C\). We note that in this regime, the main part of the energy is in the form of 
gravitational wave energy for the fast mode, whereas most of the energy is electromagnetic for the slow mode. This 
can be seen from the approximate dispersion relations above, together with the coupled equations l|39|l and (|40|l . For 
very long wavelengths, k < y/ kBq C\, we of course still have two roots of the dispersion relation, but in this case 
the two modes divide their energies roughly equal between the gravitational and electromagnetic form. Therefore we 
refer to this regime as that of mixed modes. 

The next step is to include higher order terms in the amplitude expansion of V(B X ) and let the wave amplitude 
vary weakly in space and time, as compared to u> and k. Keeping terms up to third order in B x , Eq. (|41|l now reads 

15 ^ « \ o _ kB o 



(dt+d z )[dt+[l-^ T + ^Tb~B x - jer&Bi )d z )B x = -^B x (45) 

j4 A. A 
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We note that the nonlinear terms induces second harmonic (SH) and low frequency (LF) perturbations, and we must 
therefore modify our ansatz according to: 

B x = B(z, t)e i{ - kz -^ + B SH (z, i)e 2i <**- w *> + c.c. + B LF {z, t) (46) 

where c.c. stands for complex conjugate of the preceding terms, and we define u> and k in 146|) to fulfill the dispersion 
relation l|43[> exactly. Inserting this into l|45|) we get a long equation involving both SH- and LF-terms as well as terms 
of the original frequency lo. To analyze this we use standard techniques for nonlinear wave equations, which we give 
an outline of here (see e.g. Ref. |23j for details): 

• First we note that the induced SH- and LF-terms are one order smaller in an amplitude expansion than terms 
of the original frequency, which allows us to neglect all terms of higher order in Bsh and Blf- 

• From linear theory we know that the envelope of the original perturbation B(z,t) travels with the group velocity 
v g . Since quadratic perturbations of the original amplitude acts as a driver for both second harmonic and low- 
frequency perturbations, we observe that also Bsh as well as Blf will propagate with the group velocity. Thus 
we may use the approximation dt ~ ~v g d z for derivatives acting on the amplitudes. Note also that in general 
we must use this approximation here instead of dt « — d z . This is because of the operator (dt + d z ), which 
otherwise will give zero to lowest order when acting on the wave perturbation. 

• The different time scales in the equation (second harmonics, low frequency original frequency) are picked out 
by multiplying with the appropriate exponential function and then averaging over many wavelengths and period 
times in order to omit the rapidly oscillating terms in the resultant equation, see J23|. (E.g. to find the equation 
governing the wave of original frequency, w, we multiply by exp(— i(kz — tot)) and take an average over several 
wavelengths and period times.) 



To lowest order, the equation for the SH-terms becomes 

3(ujk-k 2 ) , , x 

D(2co, 2k)B SH = \ n2 ' B 2 (47) 

where, from l|43[l . we note that D(2u>, 2k) = 3kBq/2. Similarly, the equation for the low frequency wave is 



,2 r+(A.. au, k Y\ o - 3(4a,-4fc + fc/C|) 

2B q {uj -k + k/2Cj 




where, in arriving at H48JI . we have used dt ~ —v g d z and v g is given by 

_ 2u(l-l/4C%)-2k(l-l/2C%) 



2uj - 2fc(l - 1/4C; 



(49) 



A> 



Note that if C\ — ► oo, i.e. the plasma density goes to zero, v g — > 1, although, as can be seen from l|43() . u> — k ^ 0. 
The back reaction at the original frequency is determined by the equation 

i {d t + v g d z ) + V -^dl\ B = -^r (b SH B* + B LF B + ^L-\B\ 2 b) (50) 



where the group dispersion, v' is given by 



kB 2 



(51) 



F = N (52) 



2C\ (2k-2u-* 



and the star denotes complex conjugate. 

In order to simplify the system of equations (|47|) - (|50|l . we first consider the LF-equation (|48|l . We see that for the 
different modes we can compare the solution of equation (|4*S|l with the SH-term and get Blf/ Bsh ~ d 2 \B\ 2 /k 2 . Now 
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remembering that the wave amplitudes vary weakly in space, as compared to k, we then conclude that this quantity 
is very small and we therefore neglect the term involving the low frequency held from now on. 

Let us now write equation l|50(l in a more simple form by using (|47l) to eliminate B$h and making the following 
coordinate transformations 

£ = 0(z-v g t) (53) 

W I 8 2 

r = (54) 

where 

p = C\ (4u - 4k + (55) 

and the absolute value in the transformation l|54|) has been introduced in order to keep r positive at all positive times, 
t. Note that the sign of v' g depends on the roots of the dispersion relation; for the mode with positive sign in (|44|) v' g 
is positive and for the mode with negative sign v' g is negative. Furthermore we use the linearized equations to relate 
B to the gravitational perturbation h + . Keeping only lowest order terms in Eq. I|39|) we get 



2{tu-k + ^ 



B (56) 



B uj 

Transforming the gravitational perturbation amplitude according to 

h+ = (57) 



where 



Zuj 2 k (buj - 5fc 
IQC\ (uj-k 



9k 
2CT 



n = ~ TT (58) 

k 
IT-. 



then allows us to rewrite Eq. "ST")) as a standard NLS equation for the rescaled gravitational perturbation amplitude: 



(id T ±dl)h+ = ±\h+\ 2 h (59) 

Here the plus and minus sign correspond to the fast and the slow mode respectively. Eq. (|59|l can be solved by the 
inverse scattering technique, as discussed by for example Ref. |23j . Asymptotically in time, the solution is a train of 
solitons. 



B. Generalization to a 3D spatial dependence 



In this section we address the fact that in reality we will not have exact plane wave solutions to the linearized 
Einstein and Maxwell equations. Rather the involved quantities will also depend on the x- and y-coordinates if we 
have wave propagation in the z-direction. We will here assume that this deviation from plane waves is small, and 
use a perturbative treatment. The underlying assumption is that d x and d y is of order of 1/L±, where 1/L± is the 
characteristic width of the pulse fulfilling X/L± <C 1 (A is the wavelength), while d z is of order 1/A. 

Assuming an x- and y-dependence in the metric perturbation function, will induce perturbations in the form 
of additional components in the Einstein tensor, as seen from the Lorentz gauge condition ltj(J|) . We thus note that 
there must be corresponding perturbations in the energy-momentum tensor in order to fulfill Einstein's equations 
(|19fl . The metric must of course be adjusted correspondingly, and this leads to the addition of components in the 
metric perturbation, h a b] specifically the following components (of order X/Lj_) must be added: hoi, ^02, and /i23- 
Inserting this metric perturbation into l|19|) and using the Lorentz-condition 

d b h ab = (60) 

we get the following modified version of Eq. 1(2*21 : 



(d 2 - d 2 - V\)h+ = k(8Tu - ST 22 ) (61) 
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where we again have subtracted the n and 22 components of Eq. I|19|) . and we have introduced X/\ = d 2 + d 2 . 

The variation of the perturbed fields in the x- and y-directions will also induce other components of the elec- 
tromagnetic and velocity field perturbations, so we introduce additional terms in the perturbed fields according to 
E = E x e\ + E y e2 + E z e^ and correspondingly for other fields. These extra perturbations are small, of order X/L±. 
Note that we now have several small quantities here and we neglect terms of the type d x h +1 d x E Xl v x B y etc., i.e. 
we neglect all terms that include more than one small quantity (with the important exception of the V^-operator in 
(|61H . because this operator compares with the operator (d 2 — d%), which is small in itself). We also make the same 
approximations made previously. 

With these approximations we can calculate the effective charges, currents and forces (|I3fl - l|18[) and work out the 
Maxwell and fluid equations l|7))- ljT2l . The result is a large set of equations, that can be reduced using d z « — 9 f , and 
the first order results j|57Jl and l|3*8Tl . We finally arrive at the following equation for the magnetic field perturbation 

(d 2 - d 2 - Vi) B x + ^d z ( ( Bq ^° 2B J 2 d z B x ) = B $h + (62) 

and Eq. (|oT)l reduces to 

(8 2 - 3l - Vi)ft+ = ~2 K B Q B X (63) 

In order to simplify this system further we expand the operator (d 2 — d 2 — Vj_) « 2d t (dt + d z ) — V\ and insert this 
in H62fl - I|63l) . Let us also introduce the notation <9 t _1 for the inverse operation of <9 f . Using (|63l) we can then combine 
these two equations into the following modified version of the wave equation (|41|1 

d t +d z - \dt^l) (dt + V(B x )d z - \d^V\\ B x = ~^B X (64) 

Next we note that for perturbations on the form B x ~ B(t) exp i(kz — cut) we can expand the operator d t as 
d t = —ioj + d t , where d t acts only on the amplitude. If B(t) varies slowly compared to the exponential part (i.e. 
|j- 1), we can then expand the inverse operator, 9 t _1 as 




Using this result and expanding the velocity function V(B X ) up to second order in B x we get the equation 



(dt + d z ) (8 t + V exp (B x )d z ) B x - [-■ \ d t + d z - -^rd z ) ViJ B x = ~^-B x (66) 

where V exp (B x ) = 1 - (1/2C 2 A ) + (3B X /2C^B ) - (lhB 2 x /AC 2 A Bl). We see that this is just equation ggj) with an 
extra term on the left hand side corresponding to the small variation in the x- and y-directions. In order to analyze 
this equation we make a small modification to our ansatz l|46|l in order to include the weak spatial dependence in the 
amplitudes and insert this into l|66|l . Working out the algebra, we see that the equations for the SH- and LF-terms 
are unaltered to the desired accuracy. This means that the only alteration to our previous results in section llV Al 
will be an extra term in the operator on the left hand side of equation JSUJ), corresponding to the V^-term in (|66|1 . 
Making the same transformations made in section llV Al the final result for the transformed gravitational perturbation 
amplitude, h + , is: 

(id T ± d 2 + TW\) h+ = ±\h+\ 2 h+ (67) 

where the coefficient T is given by 

( Aw - 4fc + -M- 



Auo(uj - k) (uj - fc + 2%) 



(68) 



and where the plus and minus signs again refer to the fast and slow mode respectively. 



10 



V. SUMMARY AND DISCUSSION 



We have studied the weakly nonlinear propagation of coupled Alfven and gravitational waves (AGW) propagating 
perpendicular to an external magnetic field in a plasma with C\ 3> 1, using the coupled Einstein and Maxwell 
equations. One of our main results is the general evolution equation 1|41|> . which holds for a broad band spectrum. A 
thorough study of this equation is a project for further study. In order to simplify i|41[l . we have made a WKB-ansatz, 
and showed that this leads to the well-known NLS-equation l|59|) , where the generalization to a 3D spatial dependence 
has been described in section llV Bl 

Possible astrophysical sources for large amplitude AGW are binary pulsars, particularly in the later stages of their 
evolution, i.e. not too far from merging. Clearly such a system produces gravitational waves, and, as discussed by for 
example Ref. 0], this leads also to large amplitude electromagnetic perturbations. Through the back reaction in EE 
the electromagnetic perturbation is coupled to the gravitational one, and linearly t his is described by the dispersion 
relation i|43|) . For typical wavelengths and background parameters (i.e. for k ^ y/ kBq C\), the dispersion relation 
approximately separates in the fast mode, u ~ k, and the slow mode, to « k(l — 1/2C 2 ,), where the former mainly 
has gravitational character and the later one electromagnetic character. However, very close to collapsing binaries 
we have magnetic fields B ~ 10 8 T , in which case there is no clear distinction in character between the modes. On 
the other hand, for moderate values of Bq further from the source, the modes will be clearly separated. However, 
due to the strong coupling in the near region of pulsars, both the fast and slow mode will be represented at larger 
distances, although the energy density of the later will be smaller. The amount of energy given to the slow mode can 
be estimated from the results of Ref. 0], and correspond to an electromagnetic power of the order 10 25 W close to 
merging. 

In order to discuss the nonlinear evolution, we consider pulses including a perpendicular dependence that is described 
by Eq. (|67|) . For simplicity we study long pulses with a shape that depend only on a (normalized) cylindrical radius, 
f = sj (x 2 + y 2 )/T, thus neglecting dispersive effects. On reinstating the speed of light c, the condition for this is 
\QttGB 2 C\L 2 l / c 6 /i <C k 2 L 2 , where L± and Ly are the characteristic length scales for variations in the perpendicular 
and parallel directions respectively. In this case we get a cylindrically symmetric NLS equation 

Contrary to Eq. (|59|l . there are no exact solutions known for the cylindrical version (|69fl , except for some physically 
uninteresting special cases. Still Eq. Ij69(l has been studied in some detail, both analytically, using approximate 
variational techniques |3J| , and numerically |34| . Of most interest to us is the case with the minus-sign, corresponding 
to the slow, electromagnetically dominated mode. The main motivation for considering this mode, is that this choice 
gives a nonlinearity of focusing type. For strong enough nonlinearity the solutions to l|69|) then show wave collapse, 
i.e. the pulse focuses indefinitely and the pulse radius r p — > in a finite time. The main characteristics of the collapse 
can be described within an approximate variational framework. Following Ref. |33l ] we use a trial function h + (f, r) = 
A(r)sech(f/a(r)) exp(i6(r)f 2 ) together with Rayleigh-Ritz optimization in order to derive a single ordinary differential 
equation for the pulse width a(r). Pulse energy conservation then follows from the optimization scheme according 
to a(r)A(r) = A(0)a(0). Furthermore, the evolution of the width scales as [a(r)(a(0)] 2 - 1 oc (1 - A(0) 2 a(0) 2 // c )r 2 
where I c w 1.35. Thus if A(0) 2 a(0) 2 > 1.35, the pulse will collapse to zero width in a finite time, and in the opposite 
case linear diffraction will dominate to spread out the pulse indefinitely. Of course, for the case of collapse, higher 
order nonlinearities neglected in the derivation of l|69(l will eventually become important, which will change the later 
stages of the given scenario. 

Let us study whether the electromagnetically dominated slow mode, excited by binary pulsars close to collapse may 
undergo wave collapse. Here a word of caution is at hand. In a real system of this type a number of effects outside 
the model equation (|69l) arc likely to play important roles for the pulse dynamics. For the sake of simplicity we will 
here exclude such effects. From the variational approach mentioned above, which agrees with numerical works |34| . 
we find that collapse takes place if \h + \ 2 > I c /?p, which can be written 

B 2 £ ^ok 2 c\ 

B 2 C\ SirGB 2 k r ? > (70) 

where we have reintroduced dimensional quantities. The electromagnetic fields leaving the binary system is excited 
by gravitational quadrupole radiation, which has a certain directionality, but not a very pronounced one. Thus before 
significant focusing takes place, the pulse width can be replaced by the radial distance from the source, r^istj a s a rough 
order of magnitude estimate for the pulse width. Thus we substitute r p — > rjist in i|70|) . use B 2 ^^ — 10 25 fj,o / cAtt 
Tm 2 and take k — 10~ 5 m _1 corresponding to the parameter values mentioned above, discussed in more detail in Ref. 
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[l3[ . The only extra parameters we need to specify are -Bo and po- For a broad range of relevant background values 
1(T 20 kg/m 3 < po < 10~ 5 kg/m 3 and 1(T 10 T < B < 0.1 T, the condition (|7UJ is always fulfilled. The main concern 
is to find values fulfilling the previous assumption C\/c 2 = B 2 //j,oPqc 2 3> 1, which holds for a low-density plasma 
po ~ 10~ 20 kg/m 3 , provided Bo > 1CP 4 T. However, we note that the condition l|7U|l can be misleading, since for the 
strong powers considered, the level of the second harmonic fields typically obeys Bsh > Bq at relevant distances from 
the source. This level lies outside the regime of validity for weakly nonlinear waves. We note that the left hand side 
of l|7(J|) is proportional to Bsh- Replacing the expression for Bsh found from (|47|l by its limit of validity Bsh ~ Bq 
(relevant for more moderate powers than 10 25 W), we note that collapse will occur when 



dist 



> 1 (71) 



which is easily fulfilled at reasonable distances (raist ^ 10 7 m) from the source. The characteristic timescale for 
significant focusing corresponding to the example given above (with Bsh ~ -Bo) is given by Tf oc ~ C\/c 2 uj. 

The nonlinearity of the fast (gravitationally dominated) mode is not of focusing type. We note that from Eq. Ij59(l 
there is still the possibility of other nonlinear effects, such as (dark) soliton formation that might, in principle, be 
detected. The characteristic time scale for such nonlinear processes is typically very large j^, which can be seen 
from eqs. (I54|l - (|59() inserting realistic background parameter values for the relevant quantities. This means that for 
nonlinear effects to be important for this mode, we must assume extreme parameter values which cannot be justified 
by astrophysical observations. 

The condition (|71() does not contain gravitational parameters, which is a consequence of the electromagnetic domi- 
nance of the slow mode. However, we note that the considered process is induced by the gravitational-electromagnetic 
coupling, and thus it still has a gravitational origin. The condition for wave collapse of this mode can be fulfilled for 
a reasonable range of parameters, which opens up for the possibility of structure formation of the electromagnetic 
radiation pattern. In the example considered above, the focusing takes place on fractions of a second (unless C\ is 
extremely large), and thus the nonlinearities may cause noticeable structures in the later stages of binary merging. 
We wish to point out here, though, that this example is to be seen as a somewhat crude estimation of the focusing 
effect, and that this example do not express all of the physics involved. In the absolute vicinity of a binary merger, 
other significant effects are almost certain to appear. A more complete treatment is a project for future research. 
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